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ABSTRACT 

"Green fuzzies" or "extended green objects" were discovered in the recent Spitzer GLIMPSE survey 
data. These extended sources have enhanced emission in the 4.5 /im IRAC channel images (which are 
generally assigned to be green when making 3-color RGB images from Spitzer data). Green fuzzies 
are frequently found in the vicinities of massive young stellar objects, and it has been established that 
they are in some cases associated with outflows. Nevertheless, the spectral carrier(s) of this enhanced 
emission is still uncertain. Although it has been suggested that Br a, II2, [Fe II], and/or broad CO 
emission may be contributing to and enhancing the 4.5 fiia flux from these objects, to date there have 
been no direct observations of the 4-5 /im spectra of these objects. 

We report here on the first direct spectroscopic identification of the origin of the green fuzzy emission. 
We obtained spatially resolved L and M band spectra for two green fuzzy sources using NIRI on the 
Gemini North telescope. In the case of one source, G19.88-0.53, we detect three individual knots of 
green fuzzy emission around the source. The knots exhibit a pure molecular hydrogen line emission 
spectrum, with the 4.695 fim v — 0-0 5'(9) line dominating the emission in the 4-5 /im wavelength 
range, and no detected continuum component. Our data for G19.88-0.53 prove that green fuzzy 
emission can be due primarily to emission lines of molecular hydrogen within the bandpass of the 
IRAC 4.5 /zm channel. 

However, the other target observed, G49.27-0.34, does not exhibit any line emission and appears to 
be an embedded massive young stellar object with a cometary UC HII region. We suggest that the 
effects of extinction in the 3-8 /im wavelength range and an exaggeration in the color stretch of the 
4.5 /im filter in IRAC RGB images could lead to embedded sources such as this one falsely appearing 
"green" . 

Subject headings: ISM: jets and outflows — ISM: lines and bands — ISM: molecules — ISM: 
individual(G19.88-0.53, G49.27-0.34) — shock waves — stars: formation — in- 
frared: stars — infrared: ISM 



1. INTRODUCTION 

The Spitzer Space Telescope's large scale survey of 
the Galactic plane, GLIMPSE (Galactic Legacy Infrared 
Mid-Plane Survey Extraordinaire; Benjamin et al. 2003) 
provides a rich archive of data relevant to the formation 
of stars at infrared wavelengths. This survey imaged the 
Galactic plane with the IRAC camera (Fazio et al. 2004) 
with filters at 3.6, 4.5, 5.8, and 8.0 fim at an angular 
resolution of <2" in all bands. IRAC images of mas- 
sive star formation regions often show conspicuous ar- 
eas of extended emission that appear to have enhanced 
flux in the 4.5 /xm filter (which is usually assigned to be 
the green channel when constructing 3-color images from 
Spitzer IRAC data). These regions have therefore been 
given the name "green fuzzies" (Chambers et al. 2009) 
or "extended green objects (EGOs)" (Cyganowski et al. 
2008). 

The nature of these "green" objects is not definitively 
known, but they appear to be largely associated with 
circumstellar gas and dust around, or outfiows from, 
massive young stellar objects (MYSOs). However, since 
Spitzer did not have spectroscopic capabilities covering 
the IRAC 4.5 /im filter, the specific spectral carrier(s) 
enhancing the emission within the wavelength range of 
just this one filter has not been conclusively determined. 

The IRAC 4.5 fim filter (also known as "Channel 2") 
extends from 4.0 to 5.0 ^m. If the enhanced emission 



arises from circumstellar material around the MYSO, 
lines such as Br a, Mg, and Ar, as well as the broader 
bands of deuterated PAHS at 4.48 and 4.65 //m (which 
have recently been found in ionized regions in the vicinity 
of massive stars; Peeters et al. 2004) could be the spectral 
carriers. However, if the emission arises from outflows or 
shocks, molecular and atomic shock indicators such as 
several H2 and [Fe II] lines, as well as a broad CO {v 
= 1-0) bandhead at 4.6 /im may be responsible for the 
enhanced "green" appearance of these objects. 

Several studies have tried to determine the spectral 
carriers responsible for the enhanced 4.5 /im emission via 
models or indirect observations. For instance, Noriega- 
Crespo et al. (2004) used IRS 5 — 37 /im spectra to 
study the HH47 outflow (which has particularly strong 
4.5 /im emission), and concluded based upon the H2 line 
strengths at A >5 /im that the green fuzzy emission may 
be due to H2, with some contribution by small dust par- 
ticles. Likewise, models by Smith & Rosen (2005) pre- 
dicted H2 line intensities from outflows, and showed that 
a flux enhancement would be strongest in the 4.5 /im 
fllter of IRAC. Alternatively, Marston et al. (2004) sug- 
gested that the IRAC 4.5 /im emission in the DR21 out- 
flow is due to CO emission based on the spatial coinci- 
dence of green fuzzy features to features in 115 GHz CO 
maps. 

We present here observations using Gemini Telescope's 
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near- infrared spectrometer, NIRI, of i and M band spec- 
tra of two green fuzzies from the survey of Cyganowski 
et al. (2008). The wavelength range of our observations 
cover as much of the Spitzer 3.6 and 4.5 /im fihers as can 
be observed from the ground. From these observations, 
we attempt to conclusively identify the exact spectral 
carricr(s) responsible for green fuzzy emission. 

2. OBSERVATIONS 

We selected eight northern hemisphere sources of green 
fuzzy emission with high surface brightnesses from the 
survey of Cyganowski et al. (2008); however because 
of time constraints, only the two with the highest sur- 
face brightnesses were observed: G19. 88-053 and G49.27- 
0.34. The observations were acquired with the NIRI in- 
strument on Gemini in queue mode over 7 nights in 2009. 
The f/6 camera with a 6.4 pixel slit and a plate scale 
of 0.117" pix~^ was used for all observations. The L 
and M grisms were used to obtain spectra spanning the 
wavelength ranges 2.95 — 4.12 /im and 4.47 — 5.75 ^m, 
respectively, with a spectral resolving power of ^-^460 for 
both bands. All targets were observed in "pair mode", 
in which the objects were nodded along the slit between 
two positions ("A" and "B" ). The images presented 
by Cyganowski et al. (2008) indicate that the emission 
from G19.88-053 and G49.27-0.34 have a spatial extents 
of ~30" and ~20", respectively; the separations between 
the two nod positions were therefore chosen to be 40" and 
30", respectively, to account for their extended nature 
and ensure that extended emission was not subtracted off 
when the A-B pairs were subtracted from one another. 
In both cases, the slit was intended to be set up on the 
MYSO and rotated to cut across the brightest parts of 
the extended 4.5 /im emission. The slit position angles 
were 89° in the case of G19.88-0.53, and 295° for G49.27- 
0.34. However, because neither the MYSO nor the green 
fuzzy emission knots were visible in the J band acquisi- 
tion images used to set up the observations, a blind offset 
from a star ~15" away was required to position the slit 
on the target (NIRI cannot perform spectroscopic setups 
in the M band). Inaccuracies in the absolute pointing 
and slit positioning could lead to the sources not being 
well-centered in the slit and (along with seeing fluctua- 
tions) could create significant variation in the final fluxes 
observed. Underscoring this point are the M band obser- 
vations of G19.88-0.53, which were taken over 5 nights 
and showed ^30% rms variation in the measured flux 
levels. 

All observations at L band were taken with frame times 
of 1.0 s, and the M band observations were taken with 
frame times of 0.15 s for G49. 27-0.34 and 0.2 s for G19.88- 
0.53. Coadding between 30 and 100 frames and between 
25 and 100 A-B pairs was used to reach the total integra- 
tion times. For G49.27-0.34 the total integration times 
were 1500 s at i and 1380 s at M, with all data col- 
lected on the same night. The total integration time for 
G19.88-0.53 at L was 1500 s, with all data obtained on a 
single night. At M, data were taken over 5 nights, with 
a final total exposure time of 6110 s. 

Observations of a telluric standard star and fiat field 
frames were obtained in the appropriate passbands on 
each night immediately before or after those of the target 
object. From those standard star observations the image 
quality of the observations was determined. In the case 



of G49.27-0.34, the L and M band images had a aver- 
age full width at half maximum (FWHM) of 0.69" and 
0.71", respectively. The L band image for G19.88-0.53 
had a FWHM of 0.61", and the combined data from the 
5 nights of M band observations had an effective FWHM 
of 0.51". 

The data were reduced in IRAF using the Gemini re- 
duction package specifically developed for reducing NIRI 
and GNIRS slit spectra. The processing steps included 
generating fiat field images, generating "sky" frames for 
each source frame from the observations at the opposite 
nod position, subtracting the sky frame from each image 
and dividing by the normalized fiat field, stacking the 
sky-subtracted, fiat-fielded frames and combining them. 
Residual sky emission was then subtracted from the re- 
sult. The final two-dimensional spectral images (fiux as a 
function of wavelength and position along the slit) for the 
two sources are shown in Figs. 1 and 2. The spectrum of 
the central source in each combined two-dimensional im- 
age was then traced and extracted. Multiple extraction 
apertures were defined for the various emission regions 
(including the MYSO itself) detected along the slit in 
the case of G19.88-053. For G49. 27-0.34, two apertures 
were defined corresponding to the MYSO and the weak 
extended "green" emission region. Wavelength calibra- 
tion was performed using sky emission lines. 

Telluric corrections and fiux calibration were per- 
formed in IDL using the spectra of the telluric stan- 
dards observed on each night and the general version of 
the xtellcor routine described by Vacca et al. (2003). 
For each of the various apertures defined for G19.88-0.53 
the telluric-corrected M band spectra from the sepa- 
rate observing nights were scaled to a common level (the 
scale values were derived from the observed and median 
strengths of the 4.69 /im emission line from the various 
nights) and median combined. The final L and M band 
spectra for these apertures as well as for the MYSO in 
this source are shown in Fig. 3. The final spectra for 
G49.27-0.34 are shown in Fig. 4. Line fiuxcs were mea- 
sured on the green fuzzy spectra using IDL to directly 
integrate the fiux under the features between two wave- 
length limits. 

3. DISCUSSION 

3.1. G19.88-0.53 

The extended green fuzzy emission from G19.88-0.53 
morphologically appears to be an outfiow from a MYSO 
(Fig. 1). We detect emission from three green fuzzy knots 
(GFl, GF2, and GF3) and from the MYSO itself. As 
far as we are aware, this is the first green fuzzy emis- 
sion source to be directly observed spectroscopically at 
the wavelengths covered by the Spitzer 4.5 /im channel. 
As shown in Fig. 3, we find conclusively that the main 
source of the enhanced green fuzzy emission in this source 
is the J/ = - S(9) line of Ha at 4.695 /im. We do 
not detect any other spectral features within the IRAC 
Channel 2 passband accessible from the ground; in par- 
ticular we see no emission from the putative CO feature 
at 4.6 /xm. Furthermore, we detect no continuum emis- 
sion from the green fuzzy knots. These observations of 
G19.88-0.53 demonstrate that green fuzzy emission can 
indeed be created in outfiows predominantly as a result 
of H2 line emission, without the need for CO, [Fe II], or 
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Br a emission. 

Interestingly, at the spatial resolution of our observa- 
tions (~0.7"), we detect very little spatially extended H2 
emission along the slit in our two-dimensional spectral 
images of the outflow. Any possibly spatially extended 
H2 emission must be either below our detection thresh- 
old, or located outside of the region sampled by our slit. 
This is consistent with the picture that H2 emission is 
enhanced in the tips and wakes of bow shocks formed in 
outflows (Allen & Burton 1993). 

For GFl, the position of the peak seen in the NIRI 
spectral image is spatially coincident with the peak pixel 
seen in the 4.5 /im IRAC image. GFl appears to have 
a lower surface brightness extension to the north that is 
below our sensitivity limit. The precise position of the 
GF2 peak cannot be distinguished in the IRAC image 
because of the poor IRAC resolution and the proximity 
of GF2 to the MYSO. The brightest pixel near GF3 in 
the IRAC image is actually coincident with a background 
star seen as a continuum source in our two-dimensional 
spectral image, but the extended 4.5 fim emission in the 
IRAC image does encompass the location of the knot 
detected in our NIRI spectral image. The spatial coin- 
cidence between the knots detected in our slit spectra 
and the locations of the green fuzzy emission seen in the 
IRAC images adds further support to the conclusion that 
the II2 emission we are seeing is responsible for the ma- 
jority of the enhanced 4.5 /im emission from this source. 

We attempted to extract flux densities of the individ- 
ual green fuzzy knots from the Spitzer IRAC data. How- 
ever, because of the coarse resolution of the IRAC data 
and the fact that G19.88-0.53 lies partially in a ridge of 
dust that becomes apparent only at wavelengths >5 fim, 
aperture photometry could reliably be performed only 
on the GFl knot. We used several combinations of aper- 
ture radii and sky annuli for which aperture corrections 
are available in the IRAC Handbook^. First we used 
a 2 pixel radius aperture centered on the source, with 
sky background measurement derived from the statis- 
tics in an annulus ranging from 2-6 pixels. We repeated 
these measurements with a 2 pixel aperture radius with 
a 12-20 pixel background annulus, a 3 pixel radius with 
a 3-7 pixel annulus, and a 3 pixel radius with a 12-20 
pixel annulus. The average of these four values and their 
standard deviations were calculated. 

Resultant point source flux densities for GFl are 
5.6±1.7, 18.9±6.8, and 15.8±6.3 mJy at 3.6, 4.5, and 
5.8 /im, respectively. Because the diffuse, spatially vari- 
able background dust emission is much brighter at longer 
wavlengths, the 4.5 and 5.8 /xm flux densities have a large 
uncertainties. In fact, there is no source visible at the 
location of GFl in the IRAC 8.0 /im image above this 
background. As mentioned above, GF2 is too close to the 
MYSO and GF3 lies ~2.5" from a background star, so 
reliable flux density estimates for these two sources could 
not be obtained. The background star near GF3 can be 
seen in 2MASS images at J, H and K, and demonstrates 
the need for observations with higher spatial resolution 
than afforded by Spitzer in isolating the flux from the 
green fuzzy knots themselves. 

Is the emission we are seeing from II2 lines enough to 
account for the enhancement seen in the 4.5 /im IRAC 

^ http://ssc. spitzer. caltcch.edu/irac/iracinstrunicnthandbook/32/#. 



image? As stated above, the extent of the green fuzzy 
emission seen in the IRAC images is clearly far larger 
than the width of our slit, and we are likely only detecting 
only the brighter knots of emission. With these caveats 
in mind, we converted the above IRAC photometry to 
surface brightnesses and compared them to the surface 
brightnesses derived from our NIRI spectra. The above 
IRAC flux densities translate to surface brightnesses 
of 4.8(±1.2)xl0-i'^ Wm-2/im-iarcsec-i at 3.6 /im, 
and 1.0(±0.3)xl0^^^ Wm~^/im~-^arcsec~^ at 4.5 /im. 
From the NIRI spectra we derive surface brightnesses of 
3.2(±0.2) X 10^1'^ Wm-2/im-iarcsec-i within the 3.6 /xm 
IRAC band, and 6.2(±0.4)xl0-i^ Wm-Vm-^arcsec-^ 
within the 4.5 /xm IRAC band. This latter value is de- 
rived by taking into account the estimates for the fluxes 
of the unobserved 4.180 fim ly = 0-0 S'(ll) and the 4.408 
/im f = 0-0 S{10) lines at the measured gas temperature 
and the measured ortho-to-para ratio of 1.55 (both are 
derived in §3.1.1). These data show that the II2 lines 
observed in the NIRI spectra contribute ~66% of the 
surface brightness measured in the IRAC 3.6 /xm image 
and ^62% of the surface brightness in the IRAC 4.5 /im 
image. 

We converted these surface brightnesses to flux 
densities by scaling by the size of the NIRI 
extraction aperture. This yields flux densities 

of 7.1(±2.3)xl0-i'^ Wm~2^j^-i ^^^ 35 ^ni, and 
1.2(±0.4)x 10-1*5 Wm-^m"^ at 4.5 /xm. The flux 
densities measured from the II2 lines in our spectra 
are 4.7(±0.3)xl0-i^ Wm-^/im-^ for the 3.6 fiia IRAC 
bandpass, and 7.5(±0.5)xl0^i^ Wm^^/im^^ for the 
4.5 /im IRAC bandpass. The latter value again takes 
into account the missing 4.180 and 4.408 /im lines and 
the measured ortho-to-para ratio. 

From these comparisons it is reasonable to conclude 
that the majority of the observed IRAC flux is due to 
the II2 lines seen in our spectra. Moreover, given the 
errors in both the surface brightnesses and flux densities 
derived above, we can conclude that it is possible that all 
of the green fuzzy emission is due to H2 lines to within 
our measurement errors (and this does not even take into 
account the fact that the spectral lines are the median 
flux values from across all observations which varied on 
the order of 30% from night to night; see §2). 



3.1.1. The 



conditions of the H2 emitting gas in 
G 19. 8 8-0. 5 3 



Fluxes for the H2 emission lines seen in the L and 
M band spectra of the three green fuzzies (GFl, GF2, 
and GF3) detected in our observations of G19. 88-053 
are given in Table 1. The absence of the 3.376 /im v 
= 4-3 0(3), 3.437 /im j/ = 2-1 0(5), and 3.662 /im v 
— 3-2 0(5) lines in our spectra indicates that the II2 
excitation mechanism is due to shocks and not UV flu- 
orescence. The fluorescent II2 models of Black & van 
Dishoeck (1987) indicate that these three lines should 
have strengths ^30 — 60% of that of the 3.234 /im line. 
Given the signal-to-noisc of our spectra, we can confi- 
dently rule out any substantial contributions due to flu- 
orescent excitation of H2. 

As we detect several II2 lines from GFl, we can derive a 

gas temperature from the H2 excitation diagram (see e.g., 

■"^^Mi^M^ei al. 2000). Assuming the lines are optically 
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thin, we derived column densities for each upper level 
from 

47rA/obs 



N, 
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where Jobs is the observed line intensity (W m^^ sr^^) 
and A is the radiative transition probability for each line. 
The line intensities were computed from the measured 
line fluxes and the spatial areas corresponding to the 
spectral extractions. The excitation temperature of the 
gas can be derived by plotting the values of the estimated 
column densities divided by the level degeneracies g as a 
function of the energy of the upper level E^pper (see Fig. 
5). Although we detected fewer H2 lines for the other 
two green fuzzy knots, we included the estimated col- 
umn densities from the available transitions from these 
regions in Fig. 5. As can be seen from the excitation 
diagram, the results for each of the three green fuzzies 
are consistent with a linear relation between log iVuppor/ff 
and E, 



upper • 



The slope of the line is inversely propor- 
tional to the excitation temperature, while the intercept 
yields an estimate of the total H2 column density, from 
the relation 
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where Nu^ is the total H2 column density and Z(Texcit) 
is the H2 partition function, computed using the ex- 
pressions given by Irwin (1987). From a linear fit to 
the data shown in Fig. 4 for each green fuzzy we find 
Texc^t = 2570 ± 70 K and Nh^ = 9.7(±1.0) x 10^^ 
cm-2 for GFl, Te^ctt = 2100 ± 180 K and Nh^ ^ 
4.0(±1.7) X 10^^ cm-2 for GF2, and Te^ctt = 2520 ± 80 
and Nh^ = 4.6(±0.6) x 10^^ cm'^ K for GF3. If aU three 
knots of emission are excited by the energetics of a single 
outflow from the MYSO, then the similarity of all three 
temperature values seems reasonable. We note that the 
column densities are strictly lower limits as we have not 
made any corrections for (unknown) extinction along the 
line of sight to the knots. 

We note that one line in Fig. 4, the i^ = — S(8) para- 
H2 line at 5.052 ^m, seems to be too strong compared to 
the linear, single excitation temperature relation derived 
(primarily) from the ortho-H2 lines. (The other para-Il2 
line in our spectra is the very weak v = — S(12) line at 
3.996 fim, whose line flux is highly uncertain.) Following 
Puxley et al. (2000), we used the relative strengths of the 
well-detected ortho-H2 4.695 ^m line and the 5.052 /im 
line to estimate the ortho-to-para ratio from 
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where subscript o refers to the ortho transition while p 
refers to the para transition. From the observed line 
fluxes we estimate the ortho-to-para ratio to be 1.55 ± 
0.11. 

The relative line fluxes we flnd are in reasonable agree- 
ment with those predicted by the model of Kaufman & 
Neufeld (1996) for a shock velocity of 40 km s""'^ and 
a pre-shock hydrogen density of 10^ cm~^, after taking 
into account the possible enhancement of the 5.052 /im 
line due to the low value of 7. The relative fluxes for the 
L band lines are also in reasonable agreement with those 



given by Ybarra & Lada (2009) for a shock model with 
2000 < Tgas < 4000 K. 

However, our ortho-to-para ratio value of 1.55 seems 
anomalously low for the gas temperature we have de- 
rived. Observations of II2 lines in HII54 by Neufeld et 
al. (2006) have shown that values of 1.55 are not uncom- 
mon in observations of outflows, but such low ortho-to- 
para ratios are generally found for much lower gas tem- 
peratures (-1100 K). Models by Neufeld et al. (2006) 
show that for our derived temperature of ~2600 K, we 
should find an ortho-to-para ratio of ~3. Although we 
cannot explain our low value, it should be noted that for 
IIII54, Ybarra & Lada (2009) calculate gas temperatures 
of 2000-3300 K from IRAC photometry and modeling the 
II2 line contributions within the IRAC filters. 

3.1.2. The nature of the MYSO in G19.88-0.53 

The spectrum of the MYSO powering the G19.88-0.53 
outflow is shown in Fig. 3. The presence of a deep ab- 
sorption feature attributed to solid-phase carbon monox- 
ide at 4.67 /im indicates that there is at least 20 mags 
of visual extinction toward this source (Shuping et al. 
2000). In addition to CO absorption, we detect carbonyl 
sulphide (OCS; Palumbo, Tielens, & Tokunaga 1995) in 
absorption at 4.9 /<m, as well as the "XCN" absorption 
feature at 4.62 /im (Pendleton et al. 1999). These ab- 
sorption features are seen in the spectra of very embed- 
ded MYSOs (e.g., NGC7438 IRS9; Gibb et al. 2004) and 
their presence is consistent with the lack of detectable 
flux in our spectrum at A < 3.5 /im due to high extinc- 
tion. Also present in our L band spectrum are shallow 
and wide absorption features due to methanol, centered 
at 3.53 and 3.95 /im, which have also been seen towards 
luminous MYSOs (Dartois et al. 1999). 

Although the MYSO appears to be greenish in Fig. 1, 
we do not detect any emission lines in its spectrum. All 
indications are that this is simply a deeply embedded 
source. This suggests that, in addition to being extin- 
guished at A < 3.5 /im, this source is likely to have a 
deep silicate absorption feature at 9.7 /xm (see Fig. 12 of 
Dopita et al. 2006). The presence of a deep silicate ab- 
sorption feature will depress the flux in the 8 /im IRAC 
image that was used by Cyganowski et al. (2008) as the 
red channel in their RGB composite images. A suppres- 
sion of emission due to extinction/absorption within both 
the IRAC 3.5 and 8.0 /im bands may, in part, be respon- 
sible for the greenish appearance of the MYSO itself (this 
will be further discussed in §3.3). 

Comparing the point-source photometry from the 
IRAC images to the flux levels of the NIRI spectra shows 
that we recover only approximately one sixth of the over- 
all flux from the MYSO in our spectra. The MYSO cer- 
tainly looks point-like in our spectral images, yet the 
IRAC images do appear to show much extended emis- 
sion nearby. Under the assumption that the source is 
resolved in the Spitzer IRAC images, we are able to re- 
cover the flux densities of the source in the same manner 
as for the green fuzzy, GFl; the derived flux densities 
for the MYSO are are 10.0(±l.l)xl0-i5 Wm-^^im-^ at 
3.6 /xm, and 2.1(±0.2)xl0-i^ Wm-^/xm"! at 4.5 /im in 
an area equivalent to the NIRI slit aperture. Those val- 
ues are overplotted on the spectra in Fig. 3. There are a 
few reasons that the total point source fiux is so low in 
our spectra compared to that in the IRAC images, while 



The Origin of Green Fuzzy Emission 



at the same time we recover the surface brightnesses de- 
rived from the IRAC images. First, it may be very hkely 
given the bhnd pointing described in §2 that we were not 
set up on the peak of the MYSO in the sht. Second, it 
is hkely that the source is not a point source, but has 
some extended component that constitutes a significant 
part of the measured overaU flux from the MYSO in the 
IRAC images. We suggest that a combination of these 
two influences may be at work, but are unable to quantify 
exactly how much each contributes. 

3.2. G49.27-O.34 

G49.47-0.34 is another high surface brightness green 
fuzzy. However, unlike G19.88-0.53, this source does not 
morphologically look like an outflow. Instead it is very 
similar in appearance to a cometary-shaped ultracom- 
pact HII (UC HII) region (Fig. 2). 

Although our total M band exposure time for this 
source was significantly smaller than that for G19.88- 
0.53, the time for a single nightly data set was the same 
as that for G19.88-0.53. Nevertheless, we detect no emis- 
sion from H2 or any other lines in our spectra of G49.27- 
0.34 whereas the H2 lines were readily apparent in the in- 
dividual nightly spectra of G19.88-0.53. The spectrum of 
the extended emission of the source appears to be due to 
the dust continuum from the UC HII region. Again, the 
presence of a deep CO absorption feature in the MYSO 
indicates that this source is heavily extinguished. 

Given that the emission in this source is extended 
where the putative green fuzzy is located, we can cm- 
ploy the same technique as for GFl in G19.88-0.53 
and calculate the flux densities from the IRAC data 
and compare them directly to our spectra. We de- 
rive from the IRAC data the following flux densities 
within an area equivalent to the NIRI aperture used: 
3.09(±0.02)x 10-16 Wm-^m-i for the 3.6 (xm IRAC 
bandpass, and 9.55(±0.98)xl0-i6 Wm-^/xm^^ for the 
4.5 /im IRAC bandpass. These values are plotted on the 
spectra of the green fuzzy in Fig. 5, and show that the 
continuum emission detected in our spectra represents 
all of the flux from the putative green fuzzy emission. 

A similar calculation was done for the MYSO itself 
in G49.27-0.34. Like the MYSO in G19.88-0.53, the 
derived point source fluxes from the IRAC images are 
much brighter than our spectra, indicating that we are 
likely not dealing with a point source. Therefore, we 
derived from the IRAC data the following flux densities 
from within an area equivalent to the NIRI apertures 
used: 2.15(±0.48)xl0"i^ Wm^Vm"^ for the 3.6 /xm 
IRAC bandpass, and 11.2(±3.0)xl0-i5 Wm-^m-^ for 
the 4.5 ^m IRAC bandpass. These values are plotted on 
the spectra in Fig. 5, and show that the we recover the 
IRAC flux from the extended dust continuum of both the 
MYSO and putative green fuzzy in our spectra. 

3.3. Why might compact embedded YSOs appear green? 

In the survey of Cyganowski et al. (2008), 3-color im- 
ages of the IRAC data were produced using the 3.6, 4.5 
and 8.0 fiuv data. As discussed in § 3.1.2, dust extinc- 
tion in highly embedded regions can depress the 3.6 and 
8.0 /xm emission from sources. Indebetouw et al. (2005) 
also demonstrated that the Galactic dust extinction law 



actually flattens from 5-8 y^m, rather than decreasing 
as one might expect. They also emphasized that the 
9.7 /xm absorption feature can affect the flux observed in 
the IRAC 8 /xm fiher. 

Another probable source of artificially increased 4.5 /im 
flux is discussed by Povich et al. (2007). They claim 
that many of the 4.5 /im "extended sources" can appear 
to be "green" solely due to an exaggeration in the color 
stretch used. Since the 4.5 /im filter does not encompass 
any PAH emission features, images in this filter tend to 
have much fainter diffuse emission than those in the 5.8 
or 8.0 /im bands, and to compensate for this the stretch 
is often changed to brighten the emission relative to that 
of the other bands for display purposes. 

The compact MYSO in G19.88-0.53 may appear ar- 
tificially green in the 3-color images of Cyganowski et 
al. (2008) due to a combination of these aforementioned 
effects. Since UC HII regions on average lie behind 30 
magnitudes of extinction or more (Dopita et al. 2006), it 
is possible that for these same reasons even a source like 
G49.47-0.34 may appear as an extended but artificially 
"green" source. 

4. CONCLUSIONS 

We have obtained ground-based L and Ad band spectra 
of two high surface brightness green fuzzy objects from 
the sample of Cyganowski et al. (2008). Our observations 
indicate that for the outflow source G19. 88-053, the green 
fuzzy emission is found within the outfiow itself and is 
due primarily to H2 emission lines within the IRAC 4.5 
/im band arising from discrete, spatially compact knots. 
We detect no continuum emission, or any other emission 
features. Measurements of the H2 lines in our spectra 
suggest that the knots comprise shock heated gas with 
temperatures of 2100 — 2600 K and column densities on 
the order of 4 — 10 x 10^^ cm"^. For one of the knots in 
this source, we find an unusually low H2 ortho-to-para 
ratio of 1.55. 

Our observations of G49.47-0.34 indicate that this 
source consists of a MYSO deeply embedded within an 
UC HII region. We see no sign of an outflow. Both the 
MYSO and the limb of the UC HII region appear "green" 
but exhibit only bright continuum emission; wc do not 
detect any emission lines within the IRAC 4.5 /im band- 
pass. As in the case of the MYSO in G19.88-0.53, we 
suggest that there may be reasons other than enhanced 
line emission for the "green" appearance of sources that 
are not outflows. 

Based on observations obtained at the Gemini Obser- 
vatory, which is operated by the Association of Univer- 
sities for Research in Astronomy, Inc., under a cooper- 
ative agreement with the NSF on behalf of the Gemini 
partnership: the National Science Foundation (United 
States), the Science and Technology Facilities Coun- 
cil (United Kingdom), the National Research Council 
(Canada), CONICYT (Chile), the Australian Research 
Council (Australia), Ministrio da Cincia e Tecnologia 
(Brazil) and Ministerio de Ciencia, Tecnologa e Innovacin 
Productiva (Argentina). Gemini program idcntiflcation 
number associated with the data is GN-2009B-Q-79. We 
would like to thank Andrew Stephens of Gemini Observa- 
tory for his help with aspects of the NIRI data reduction. 

Facilities: Gemini:Gillett (NIRI) 
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The Origin of Green Fuzzy Emission 



TABLE 1 
Observed H2 Line Fluxes for G19.88-0.53 Green Fuzzies 
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GF2 
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Transition 


Line Flux 


Error 


Line Flux 


Error 


Line Flux 


Error 




(mid) 


(W m-2) 


(W m-2) 


(Wm-2) 


(Wm-2) 


(W m-2) 


(W m-2) 


u=l-0 0(5) 


3.234 


1.92E-17 


1.02E-18 






9.20E-18 


5.28E-19 


,y=l-0 0{7) 


3.807 


8.32E-18 


l.lOE-18 










,y=0-0 5(13) 


3.847 


1.12E-17 


1.24E-18 


1.54E-18 


5.14E-19 


4.74E-18 


5.67E-19 


,y=0-0 5(12) 


3.996 


4.97E-18 


9.87E-19 










^=0-0 5(9) 


4.695 


3.13E-17 


1.23E-18 


7.42E-18 


4.87E-19 


1.02E-17 


6.16E-19 


u=0-0 5(8) 


5.052 


2.10E-17 


1.18E-18 


2.86E-18 


3.79E-19 


3.44E-18 


5.01E-19 



Note. — The 4.695 and 5.052 fim line fluxes were measured on a median spectrum generated from 
six M band data sets taken over 5 nights. The L band measurements are from a single data set taken on 
one night. Errors are purely statistical. Absolute errors are expected to be on the order of ~30%. 




H 2 1-0 0(5) 



De Buizer and Vacca 

H20-0S[13] 
H2 1-0017)1 H2 0-0S[121 



CO 



HjO-OSfQl H 2 0-05(8) 



10" 




Green Fuzzv 2 



-10" 



3.2 



^m 



3.4 



3.6 



3S 



4.0 



Al 



AA 



i& 



5.0 urr 



Spitzer 3,6 micron channel 



jpitzer^,^' micr'jn channel 



Fig. 1. — Green fuzzy emission from G19.88-0.53. (Left) A 3-coIor composite image from Cyganowski et al. (2008), with 3.6 fj,m as blue, 
4.5 fim as green, and 8.0 fim as red. The area containing extended green emission is encompassed in the green contour. The NIRI slit 
position and angle used for spectroscopy is shown by the white line. (Right) 2-D images of the combined L and M band spectra (with 3-pixel 
smoothing for display purposes only) are shown with the same vertical spatial scaling as the figure on the left. White lines encompass the 
spectra of the three green fuzzy knots detected. The bandpasses of the IRAC 3.6 and 4.5 /im filters are shown at the bottom. The only 
detectable emission from the three green fuzzies are lines of shock excited H2, which are labeled at the top. The absorption feature of CO 
is also marked. No data are present from 4.10 to 4.55 fj,m where the atmosphere is opaque. 




Fig. 2. — Green fuzzy emission from G49.27-0.34. (Left) A 3-color composite image from Cyganowski et al. (2008), with 3.6 fj,m as blue, 
4.5 fim as green, and 8.0 fiai as red. The area containing extended green emission is encompassed in the green contour. The NIRI slit 
position and angle used for spectroscopy is shown by the white line. (Right) 2-D images of the combined L and M band spectra with 
3-pixel smoothing for display purposes only) are shown with the same vertical spatial scaling as the figure on the left. The only detectable 
emission from the extended "green" emission area is continuum. The absorption feature of CO is marked, as well as a region that results 
from poor sky subtraction (encircled cross symbol). No data are present from 4.10 to 4.55 fim where the atmosphere is opaque. 
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Fig. 3. — Combined L and M band spectra extracted from the locations within G19.88-0.53 marked as circles in Figure 1. The light 
green area delineates the passband of the IRAC 4.5 /^m filter. The only detectable emission features from the three green fuzzies are lines 
of shock excited H2, which are delineated by vertical dashed blue lines with the transitions labeled at the top. Areas of the spectra with 
missing data were removed because of poor sky subtraction. In some cases bad pixels create false features which are delineated by the 
encircled cross symbol. Flux densities derived from the IRAC image photometry of the MYSO are plotted as red circles. No data are 
present from 4.10 to 4.55 ^m where the atmosphere is opaque. 
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Fig. 4. — Combined L and M band spectra for G49.27-0.34. The light green area delineates the passband of the IRAC 4.5 fim filter. 
The only detectable emission features from the putative green fuzzy area is continuum emission. The expected lines of shock excited H2 , 
which are delineated by vertical dashed blue lines with the transitions labeled at the top. Areas of the spectra with missing data were 
removed because of poor sky subtraction. In some cases bad pixels create false features which are delineated by the encircled cross symbol. 
Flux densities derived from the IRAC image photometry of the MYSO and the putative green fuzzy are plotted as red circles. No data are 
present from 4.10 to 4.55 /im where the atmosphere is opaque. 
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Fig. 5. — Excitation diagram of the H2 level column density distribution towards the green fuzzy knots in G19.88-0.53. We are able to 
fit these data with single temperatures value for each of the green fuzzy knots (solid lines). 



